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Abstract We examine the effect on cosmological evolution of adding a 
Gauss-Bonnet term to the standard Einstcin-Hilbcrt action for a (1 + 3) + d 
dimensional Friedman-Robertson- Walker (FRW) metric. By assuming that 
the additional dimensions compactify as a power law as the usual 3 spa- 
tial dimensions expand, we solve the resulting dynamical equations and find 
that the solution may be of either de Sitter or Kasner form depending upon 
whether the Gauss-Bonnet term or the Einstein term dominates. 
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1 Introduction 
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, Ever since the introduction of extra dimensions by Kaluza and Klein, the 

notion that hidden extra dimensions may play a role in the dynamics of our 
usual four dimensional spacetime has received considerable attention (for a 
. review see Applequist et. al. [T]), [2]. However if these extra dimensions exist, 
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experiments do provide constraints on the maximum current size (of extra 
dimensions) to be ^ 100 /Ltm. [3] An idea that has been extensively studied 
is that these extra dimensions were once large but underwent dynamical 
compactifaction as the usual three spatial dimensions grew. This has been 
studied by Paul and Murkherjee |4j and Mohammedi [5] among others. 

It is also widely thought that Einstein gravity is only a low energy ef- 
fective field theory which requires modification at higher energy. A detailed 
review of this idea was given by Deser in [6]. One possible modification of 
the Einstein-Hilbert action comes from adding additional Lovelock terms. 
[7] This modification is attractive because it yields second-order, divergence 
free field equations as one would demand from a generally covariant theory of 
gravitation. 'Si One may formulate Lovelock gravity theory as an expansion 
in powers of the curvature to obtain a zeroth-order constant term or cos- 
mological constant, a first-order term that gives the usual scalar curvature 
which yields Einstein gravity, and a second-order term that is known as the 
Gauss-Bonnet term plus higher order terms. 

This paper incorporates both of these ideas: compactification of the higher 
dimensions and the addition of a Gauss-Bonnet term to the action. We con- 
sider a dynamical compactification of a D-dimensional manifold to a max- 
imally symmetric manifold of dimension d and an expanding FRW space- 
time of dimension 4 where we have modified the Einstein-Hilbert action by 
including a Gauss-Bonnet (GB) term. [9] This Gauss-Bonnet term can be 
interpreted as being a first order correction from string theory or simply a 
modification of Einstein gravity. 

We will choose a power relation between the conformal factors of the three 
spatial dimensions and the d extra dimensions. This choice is motivated by 
the fact when the model is void of Gauss-Bonnet terms, the resulting field 
equations become exactly that of 4D FRW cosmology for arbitrary values of 
n and d, after one redefines the coupling and cosmological constant. [5] 

This paper is organized as follows. In section [2] we present a general 
Einstein-Hilbert plus Gauss-Bonnet action in _D-dimensions where the field 
equations are calculated and the correction to the FRW equation is given. 
In section [3l we define an "effective" pressure which an observer constrained 
to only the four dimensional volume would observe. This arises from the 
conservation equation, where this defined effective pressure yields an identical 
expression for the conservation equation to that of a 4D theory. In section 
m we find specific solutions for the A + d dimensional Einstein plus Gauss- 
Bonnet equations for cases when the Einstein terms or the Gauss-Bonnet 
terms dominate. Finally in section [5] we summarize our results. 



2 Framework for Field Equations 

In this paper we will follow the notation of Paul and Mukherjee[3] and Mo- 
hammedi [5] to express the Einstein-Hilbert action with an extra Gauss- 
Bonnet term as 




(1) 
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where ^ is a Gauss-Bonnet term {G = RabcdR-^^^^ ~ ARabR"^^ +R'^). A 
variation of the action ([T]) with respect to gAB produces the equation; 

+ A5A^ + e = -Ta^ (2) 

K 

where k is the 4 + d dimensional couphng constant and the Einstein and 
Gauss-Bonnet tensors, respectively, are 

Ga^ = Ra^ - ^A^R (3) 

Ga"" = i - 4i?BI5i?^^ + R') Sa"" 

^{^RbdeaR^^^^ + 2RRa^ — 4:Rd^ Rba^'"' — "^Rb^ Ra^) -(4) 

Furthermore we will assume that the stress-energy tensor will be that of 
a perfect fluid, thus it is of the form 

T^"" = dia.g[p{t),p{t),p{t),p{t),pd{t),...pd{t)] (5) 

where Pdit) is the pressure on the higher dimensional compact manifold. 
We will choose a metric ansatz: 



ds^ = -dt^ + a^{t) 



{dO^ + sin^ Odcf") 



J- Kr^ 

+ &'(i)7™(y)dy™dj/" (6) 

where the extra dimensions are defined to be maximally symmetric such 
that the Riemann tensor for 7 has the form Rabcd = k{'yac^bd — ^ad'Jbc) ■ 
In agreement with current observations we will consider the usual 3 spatial 
dimensions to be flat {K = 0) and also demand that the extra dimensions 
be flat (fc = 0) as well in agreement with Mohammedi [5] that one finds 
unphysical properties for p and p if fc 7^ . 

The metric leads to Riemann tensors of the following form (where both 
the dimensions are flat) 

RaOaO = I Rabab = 0,^0^ , RmOmO = bb , 
Rmama = aabb , Rabab = (7) 

where a, b are indices which run from 1...3 and to, n are indices which are in 
the extra dimensions. The Ricci Tensor and Ricci Scalar are 



a 



.b „ .. ,aab 



RoQ = 3 h , Raa — 2a + a'd + d , 

a b b 

Rm,n^3^ + id-l)b^+bb (8) 



a b ab a? , ,r 
\- + 2d- + M—+% 
a b ab 

where d is the number of extra dimensions 



R^Q- + 2d-- + &d—+&^ + d{d^l)— (9) 
a b ab a'^ b^ 
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The full Einstein tensor may be expressed as having an Einstein term 
G^'^ and a Gauss-Bonnet term G^i'^ ■ We make the assumption that the extra 
dimensions compactify as the 3 spatial dimensions expand as was done by 
Mohammedi [S] 

bit) (10) 

where n > in order to insure that the scale factor of the compact manifold 
is both dynamical and compactifies as a function of time. With this ansatz, 
the non-zero elements of the Einstein tensor take the form 



Gm 
m — — 

a n 



(2771 + 37^2) - + [2 (2771 + 37^2) - dn (771 + 3772)] -j 



(11) 
(12) 

(13) 



where we defined the coefficients 



771 = — [6 + dn{dn — n - 

772 = (d 71 - 2) 



6)] 



(14) 



Note that when D = 4 (or d = 0) the equations become the well known FRW 
equations in four dimensions. In 4 -I- c? dimensions, the Gauss-Bonnet terms 
become 



r ^ « 



^ r, 1 , > N a ^ ad 

Ga" = T7 (?i - dn 6) ^ + 6 ^ 



n 1 (in 

g^n"^ = - (6 + 6) ^ + (3 6 + 4 6) ^ 

dn a'^ 



(15) 
(16) 
(17) 



where we defined the coefficients 

^1 



6 

6 



-(i7i 



(d- 1)71 



(d - 2)77 [(d - 3)77 - 12] + 18 
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-2di 



{d - 1)77^ [{d - 2)77 - 6] + 6 



(18) 



where the constants 77^ and depend upon the values of 77 and d as defined 
above. Note that if we force the number of extra dimensions to zero (d — > 0) 
then the Gauss-Bonnet terms vanish for fJo° and Qa"^ as one would expect in 
four dimensions. 
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3 Effective Pressure and the Field Equations 

For the case of both spaces having flat curvature, the D-dimensional Fricdmann- 
Robertson- Walker (FRW) equations ^ take the form 



2k 
P_ 

P£ 

2k 



0? 

m— 

a 



/ 2\ 



1 

d n 



(2771 + 3772) - + [2 (2771 + 3772) -dri (771 + 3772)] ^ 



, a 1 , , aa 

6 + 6 4 ^1 + 3^3 — 

an a'^ 



(19) 
(20) 

(21) 



where we have set A = 0. Together with these Einstein equations, we also 
demand that the conservation equation hold for the stress-energy tensor 
(V^T^, = 0) or 



d 

At' 



Using the assumption that b ~ 1/a", this becomes 

which by simple algebra may be written in the more familiar form 

p + 3-(p + 75) = 0. 
a 



(22) 



(23) 



(24) 



Note that this is simply a statement that dE = —PdV where we have 
defined an "effective" pressure (p) [5] which an observer constrained to exist 
only upon the "usual" 3 spatial dimensions would see as 



p = p- -dn [p+Pd]- 



(25) 



As was pointed out by Mohammedi [5], this effective pressure can be nega- 
tive for positive values of p, p, and pd- The effective pressure can be easily 
computed from the d-dimensional FRW equations P^ - (PT|) and is given by 
the relation 



2k 



1 



(26) 



Note that the above field equations (fT9)) - ((2T|l . ((26| become the same as j5] in 
the limit where e — *■ (i.e. no Gauss-Bonnet terms). By redefining the cou- 
pling and cosmological constant in eqns. (fT9|) and ([26|) . one recovers standard 
4-D FRW cosmology for arbitrary values of n and d in the e ^ limit. 
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4 The General Solutions of the GB FRW Field Equations 

The effective Z?-dimensional FRW equations and the conservation equation 
now read 






= p + 3 


p_ 


= m— 


P_ 

2k 


i 


Pd_ 


f 


2k 


d n 



^ = --r;i 2- + - --66 4--- - 



(2771 + 3772) - + 2^ -dn (771 +3772)^ 



a 



d n 



3 (a + 6) ( - - 



6^ 

a 



(27) 
(28) 
(29) 

(30) 



One can easily show that ([27|l is trivially satisfied when (p8|l and ([29|) are 
employed. This is exactly like that of 4D FRW cosmology where the variables 
are left underdetermined by the field equations and one usually chooses an 
equation of state (EoS) of the form 



p = w p 



(31) 



to proceed. Notice that the parameter w can in general be time-dependent. 
If w is in fact constant, (P7|) can be integrated yielding 



p{t) = poa-3(i+'") 



(32) 



Using (|3ip and eliminating p and p from the field equations, we obtain 
an expression of the form 



(33) 



The behavior of the scale factor is determined by this expression. As is ob- 
vious from ()33|) . the scale factor is dependent on the value of the parameter 
w as one would expect. In the following subsections, we will examine the 
solutions of this equation for cases when parameter is uniquely w — —1 and 
the general case, when w ^ ~1. 



4.1 Solutions for w — —1 

If Tjj = — 1, the equation of state ([55)1 reduces to 



(34) 
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This equation has two solutions depending on which bracket is equated to 
zero. In either case, we find a de Sitter-type solution given by the expressions 



a(t) = oq e 



(35) 



where H is the Hubble constant which is obtained from the remaining field 
equations. Plugging a{t) into (|28|) yields a value for the Hubble parameter of 
the form 

2 



H 




(36) 



Notice that when the vacuum-energy density is chosen such that the radical 
vanishes, we find a result corresponding to what one obtains when the first 
bracket is set equal to zero. Hence, the first solution is simply a special case 
of the second solution. 

In the small epsilon limit, when the Gauss-Bonnet contribution offers a 
small correction to the Einstein-Hilbert term, we find 



P 



2Krji 



1 ^1 P 



(37) 



to 0{e). Note that we kept only the negative root of (|36p so that the value 
of the Hubble constant approaches that of 4D FRW cosmology in this small 
e limit. 

In the large epsilon limit, when the Gauss-Bonnet term dominates over 
the Einstein term, we find a value for the Hubble constant of the form 



±- 



1 



(3 V 2kt]i 



1 

2(3^ 



where we defined the parameter 



/3 = 




(38) 



(39) 



and kept terms to ©(l/e). 

We can eliminate the expansion factor and obtain a higher dimensional 
equation relating both pd and p. This EoS-like expression takes the form 



{Pd - ^ap) + 77 



2 i^a - xo) 



(40) 



where we defined the constants 



1 



^0 = 

m 
1 

Ti 



1- 

an: 



Xo 



1 



(a +6) 



3?72) - m 
1 



(41) 



Notice that (j40p is inherently non-linear as one might expect from an exact 
solution to Gauss-Bonnet FRW cosmology. Also notice that when dn = 3, 
which corresponds to a constant volume d + 4 dimensional Universe, (I4ip 
dramatically simplifies yielding the fixed parameters ipo = I and xo — 1/3. 
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4.2 Solution for w 7^ -1 



In this subsection, we investigate the behavior of the scale factor for the 
general case when it; is left as a free parameter. By integrating the equation 
of state, eqn. (l33l) 



a 



tan 



-(—--{1 

13 \a/a 2^ 



w)t 



(42) 



where we again defined /3 = \f^^^Jr\\■ In this section, we can explore this 
behavior by investigating (|42)) when the Gauss Bonnet term is small and 
when it dominates over the Einsteinian term. 



^.2.1 Small epsilon regime 

In the small e regime, one can expand the equation of state (|42|) by taking 
the arctan of each side and taking the expansion for small x. Keeping only 
the lowest order contribution in the expansion, it becomes 

0. (43) 



a" 3 , s a, 
(3' — + -{l + w)t-~l 
a'^ Z a 



Solving for the Hubble parameter, we obtain 
d 3 



4/3 



1 ± 



4/3 



3(1 + 'w)t 



(44) 



Again, expanding the square root in the small e 0^) limit and keeping the 
first order e contribution, one obtains 



Hit) = 



1 



a 3(l + w)t 



1 



6 



1 



3(1 + w) t 



(45) 



where we again kept only the negative root in order to arrive at 4D FRW 
cosmology in this small e limit. Using this result, we integrate (|45p and obtain 
a value for the scale factor of the form 



i{t) ~ ^i 1 2/3(1+-) 



1 



2771 



3(1 + u>) 



(46) 



to first order in e. Note that in the large time limit (t 00) the scale 
factor tends to its zeroth-order value. Hence, we see that the Gauss-Bonnet 
contribution becomes vanishingly insignificant for late cosmological times as 
one would expect. 

As a consistency check, we insert (15^ . (HSI) . and into (|28p and find 
that the time dependence does in fact vanish to first order in e. We obtain a 
relation on the coefficients given by 

l/3(l+ii)) 



-il + w) 



2^7/1 



(47) 



which agrees with the value one obtains from 4D FRW cosmology. 
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Large e regime 



Interestingly, because the first order term in the the expansions of tan(2:) and 
arctan(x) are both hnear in x\ the expansion for when the Gauss-Bonnet term 
dominates (large e or equivalently when t <ti /3). (|42p is the same as for small 



a p \a/a 2 



w)t 



(48) 



The difference between the large and small e cases is that in the large e case 
one does not expand the radical as was done in (|45p , instead one may simply 
express the equation in terms of H as 



it , 

2^(1 + -) 



(49) 



If solves this quadratic and keeps the terms large in 1/e (the terms that 
dominate when GB is large), one obtains 



a{t) « aoExp 



it 



(3 8/3' 



13^ 



(50) 



Note that when the Gauss Bonnet terms dominate the form of a{t) is de 
Sitter like. 



5 Conclusion 

In conclusion, we have studied the Friedman equations for an Einstein plus 
Gauss-Bonnet action in 4 -I- d dimensions. We furthermore demanded that 
the extra dimensions compactify as a{t) ^ 6(i)~" where n > 0. In Section^ 
we solved the field equations for cases when w = — \ and w ^ —1. We found 
solutions for both of these cases when the Gauss-Bonnet term dominates 
the Einstein term early time cosmology and when the Gauss-Bonnet term is 
subdominant to the Einstein term (late time cosmology) . We have found that 
when Einstein gravity dominates, the conformal factor is a Kasner solution 
with a small Gauss-Bonnet correction. However, when the Gauss-Bonnet 
term dominates, a de Sitter type solution is obtained indicating that the 
Gauss-Bonnet term gives rise to nontrivial corrections of the scale factor 
a{t). One may conclude from this analysis that if the Gauss-Bonnet term is 
dominant, the compactifying extra dimensions can be thought of as playing 
the role of a "cosmological constant" forcing a{t) to behave as a de Sitter 
solution even though no such constant is present in the action. 

The natural question which arises from this paper is when does this model 
of compactifing internal dimensions with extra Gauss-Bonnet terms physi- 
cally representative of our universe. Without further modification of the pa- 
per, it seems to be applicable only to early universe inflation as our model 
tends towards a Kasner solution in late time. In order to agree with current 
cosmological observations, the model should allow the extra dimensions to 
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compactify much slower then the three spatial dimensions grow to be in order 
to agree with the current limits on the running of the coupling constants. 

This paper leads to several intriguing questions. First, does this same 
effect of compactifing extra dimensions leading to an effective cosmological 
constant exist in alternate models such as brane models like that of Randall- 
Sundrum [TT] or Narain et al ,22J. If so, then could this be used to generate 
an inflationary cosmology without the use of a scalar field [13] . A second and 
central issue to explore would be to relax the explicit relation between a(t) 
and b{t), could one then solve the resulting field equations? A possibility we 
are examining seems to require a numerical study. 
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